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We discuss the perturbative QCD description of diffractive deep-inelastic scatter- 
ing, and extract diffractive parton distributions from recent HERA data. The 
asymptotic colHnear factorisation theorem has important modifications in the sub- 
asymptotic HERA regime. In addition to the usual resolved Pomeron contribution, 
the direct interaction of the Pomeron must also be accounted for. The diffractive 
parton distributions are shown to satisfy an inhomogeneous evolution equation, 
analogous to the parton distributions of the photon. 



1. Introduction 

A notable feature of deep-inelastic scattering is the existence of diffractive 
events, "f*p X +p, in which the slightly deflected proton and the cluster 
X of outgoing hadrons are well-separated in rapidity. At high energies, the 
large rapidity gap is believed to be associated with 'Pomeron', or vacuum 
quantum number, exchange. The diffractive events make up an appreciable 
fraction of all (inclusive) deep-inelastic events, 7*p X. We will refer to 
the diffractive and inclusive processes as DDIS and DIS respectively. The 
recent improvement in the precision of the DDIS data [1~3] allow improved 
analyses to be performed and more reliable diffractive parton density func- 
tions (DPDFs) to be extracted. In this article we criticise the conventional 
extraction of DPDFs based on 'Regge factorisation' in which the exchanged 
proton is treated as a hadron-like object. We show using perturbative QCD 
(pQCD) that the treatment of diffractive PDFs has more in common with 
the photon PDFs than with the proton PDFs. 



*To appear in the proceedings of the Ringberg Workshop on "New Trends in HERA 
Physics 2005", Ringberg Castle, Tegernsee, Germany, 2—7 October 2005. 
tfn collaboration with A.D. Martin and M.G. Ryskin. 
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2. DifFractive parton distributions from Regge factorisation 

Let the momenta of the incoming proton, the outgoing proton, and the pho- 
ton be labeUed p, p', and q respectively; see Fig. 1(a). Then the basic kine- 
matic variables in DDIS are the photon virtuality, — , the Bjorken-x 
variable, Xg = /{"^P • l), the squared momentum transfer, t — (p — p')^, 
the fraction of the proton's light-cone momentum transferred through the 
rapidity gap, xp = 1 —p'~^/p'^, and the fraction of the Pomeron's light-cone 
momentum carried by the struck quark, /? = Xs/xf. 




Figure 1. (a) Resolved Pomeron contribution in the 'Regge factorisation' approach, (b) 
Resolved Pomeron contribution in the 'perturbative QCD' approach, (c) Direct Pomeron 
contribution in the 'perturbative QCD' approach. 



It is conventional to extract DPDFs from DDIS data using two levels 
of factorisation. Firstly, collinear factorisation means that the diffractive 
structure function can be written as [4] 



F^'-'\xr,f3,Q^)= C2,a®a^, (1) 

o.=q,g 



where the DPDFs a° = or zg° , with z e [/3, 1], satisfy DGLAP evolu- 
tion: 

and where 6*2,0 and Paa' are the same hard-scattering coefficients and split- 
ting functions as in inclusive DIS. The factorisation theorem (1) applies 
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when Q is made large, therefore it is c;orroct up to power-suppressed cor- 
rections. It says nothing about the mechanism for diffraction, which is 
assumed to reside entirely in the input DPDFs fitted to data at a starting 
scale Qq] see Fig. 1(a). 

In a second stage [5] Regge factorisation is usually assumed, such that 



where the Pomcron PDFs a' = z(f or zg^ . The Pomeron flux factor /p is 
taken from Regge phenomenology. 



Here, ar{t) = ap{0) + a'pt, and the parameters Br, ap(0), and a'p should 
be taken from fits to soft hadron data. Although the first fits to use this 
approach assumed a 'soft' Pomeron, ap(0) 2± 1.08 [6], all recent fits require 
a substantially higher value to describe the data. In addition, a secondary 
Reggeon contribution is needed to describe the data for xr > 0.01. This ap- 
proach is illustrated in Fig. 1(a), where the virtualitics of the t-channcl par- 
tons are strongly ordered as required by DGLAP evolution. The Pomeron 
PDFs aF are parameterised at some arbitrary low scale Qq, then evolved 
up to the factorisation scale, usually taken to be the photon virtuality Q^. 

Although this approach has been found to give a good description of 
the DDIS data [1,2,7,8], it has little theoretical justification. The 'Regge 
factorisation' of (3) is merely a simple way of parameterising the xp depen- 
dence of the DPDFs. Note, however, that the effective Pomeron intercept 
ar{0) has been observed to depend on [3], contrary to the 'Regge fac- 
torisation' of (3). The fact that the required ap(0) is greater than the 
'soft' value indicates that there is a significant perturbative QCD (pQCD) 
contribution to DDIS. 

3. DifFractive parton distributions from perturbative QCD 

In pQCD, Pomeron exchange can be described by two-gluon exchange, two 
gluons being the minimum number needed to reproduce the quantum num- 
bers of the vacuum. Two-gluon exchange calculations arc the basis for the 
colour dipole model description of DDIS, in which the photon dissociates 
into qq or qqg final states. Such calculations have successfully been used to 
describe HERA data. The crucial question, therefore, is how to reconcile 
two-gluon exchange with coUincar factorisation as given by (1) and (2). Are 
these two approaches compatible? 



a {xp,z 



Q^) = fr{xr)a^{z,Q^) 



(3) 




(4) 
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Generalising the qq or qqg final states to an arbitrary number of par- 
ton emissions from the photon dissociation, and replacing two-gluon ex- 
change by exchange of a parton ladder, we have diagrams like that shown 
in Fig. 1(b) [9-12]. Again, the virtualities of the t-channel partons are 
strongly ordered: /iQ -C • . • <C /i^ <C . • • ^ Q^. The scale /i^ at which 
the Pomeron-to-parton splitting occurs can vary between /Xq ~ 1 GeV^ and 
the factorisation scale Q^. Therefore, to calculate the inclusive diffractive 
structure function, -F^^^\ we need to integrate over /x^: 



Ju.1 M 



Here, the perturbative Pomeron flux factor can be shown to be [12] 

1 



,2 

/p(xp;/i^) = 



xpB 



D 



Rg xrg{xr,fj, ) 

A* 



(6) 



The diffractive slope parameter Bd comes from the f-integration, while the 
factor Rg accounts for the skewedness of the proton gluon distribution [13]. 
There is a similar contribution from sea quarks, where g{xp,iJ,^) in (6) is 

replaced by S'(.Tp,/i^), together with an interference term. In the fits pre- 
sented here, we use the MRST2001 NLO gluon and sea-quark distributions 
of the proton [14]. The Pomeron structure function in (5), -F2 (/?, Q^; i-i^), is 
calculated from Pomeron PDFs, a'(z, Q^; ^^), evolved using NLO DGLAP 
from a starting scale /i^ up to Q^, taking the input distributions to be 
LO Pomeron-to-parton splitting functions, {z , ji^ ; ji'^) = Pav{z) [11,12]. 
At first glance, it would appear that the perturbative Pomeron flux factor 
(6) behaves as /p(.xp;/i^) ^ 1/m^i so that contributions from large ji^ arc 
strongly suppressed. However, at large /U^, the gluon distribution of the pro- 
ton behaves as xrg{xp, fx^) ~ (m^)'"') where 7 is the anomalous dimension. 
In the BFKL hmit of xp ^ 0, 7 ~ 0.5, so /p(a;p; y?) would be approxi- 
mately independent of /i^ . The HERA domain is in an intermediate region: 
7 is not small, but is less than 0.5. In Fig. 2(a) we plot (6) multiplied by jj? 
to show that (6) does not behave as l//i^ at small xr- It is also interesting 
to plot the integrand of (5) as a function of /i^, as shown in Fig. 2(b). 
Notice that there is a large contribution from /^^ > 2-3 GeV^ , which is the 
value of the input scale Qg typically used in the 'Regge factorisation' fits 
of Sect. 2. Recall that fits using 'Regge factorisation' include contributions 
from /i^ < Qq in the input distributions, but neglect all contributions from 
jj? > Qq\ from Fig. 2(b) this is clearly an unreasonable assumption. 
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As well as the resolved Pomeron contribution of Fig. 1(b), we must also 
account for the direct interaction of the Pomeron in the hard subprocess, 
Fig. 1(c), where there is no DGLAP evolution in the upper part of the 
diagram. Therefore, the diffractive structure function can be written as 

F^^'^ = J2 ^2,a a° + C2^r ; (7) 

a—q^g ^"^"^ 

V y Direct Pomeron 

Resolved Pomeron 

cf. (1) where there is no direct Pomeron contribution. The direct Pomeron 
term, C2,p, calculated from Fig. 1(c), will again depend on fp{xp; /i^) given 
by (6). Therefore, it is formally suppressed by a factor l//i^, but in practice 
does not behave as such; see Fig. 2(a). 

The contribution to the DPDFs from scales /i > /iq is 

/■"^'d/z^ 

aD(xp, z, Q^)= ^ Mxr, m') a^{z, A^')- (8) 

Differentiating (8), we see that the evolution equations for the DPDFs 
are [12] 

daP 

Paa'®a'^ + Par(z) h{^r-Q''); (9) 

o.'=q,g 

cf. (2) where the second term of (9) is absent. That is, the DPDFs satisfy an 
inhomogeneous evolution equation [10,12], with the extra inhomogeneous 
term in (9) leading to more rapid evolution than in the 'Regge factorisation' 
fits described in Sect. 2. Note that the inhomogeneous term will change the 
xp dependence evolving upwards in Q^, in accordance with the data, and 
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unlike the 'Regge factorisation' assumption (3). Again, the inhomogeneous 
term in (9) is formaUy suppressed by a factor but in practice does 

not behave as such; see Fig. 2(a). 

Therefore, the diffractive structure function is analogous to the photon 
structure function, where there are both resolved and direct components 
and the photon PDFs satisfy an inhomogeneous evolution equation, where 
at LO the inhomogeneous term accounts for the splitting of the point- 
like photon into a qq pair. If we consider, for example, diffractive dijet 
photoproduction, there are four classes of contributions; see Fig. 3. The 
relative importance of each contribution will depend on the values of a;-y, 
the fraction of the photon's momentum carried by the parton entering the 
hard subprocess, and zp, the fraction of the Pomeron's momentum carried 
by the parton entering the hard subprocess. 



Resolved photon 

{X-y < 1) 



Direct photon 

(•^7 = 1) 



Resolved 
Pomeron 
{zr < 1) 





Direct 
Pomeron 

{zr = 1) 





Figure 3. The four classes of contributions to diffractive dijet photoproduction at LO. 
Both the photon and the Pomeron can be either 'resolved' or 'direct'. 
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4. Description of DDIS data 

A NLO analysis of DDIS data is not yet possible. The direct Pomeron 
terms, C2,p, and Pomeron- to-parton splitting functions, Pav, need to be 
calculated at NLO within a given factorisation scheme (for example, MS). 
Here, we perform a simplified analysis where the usual coefficient functions 
C2,a and splitting functions Paa' {a, a' = q,g) are taken at NLO, but C2,p 
and PaP are taken at LO [12]. We work in the fixed flavour number scheme, 
where there is no charm DPDF. Charm quarks are produced via 7*5' — > cc 
at NLO [15] and 7*P — > cc at LO [16]. For Hght quarks, wc include the 
direct Pomeron process j^f qq at LO [12], which is higher-twist and 
known to be important at large (3. 

To see the effect of the direct Pomeron contribution and the inhomoge- 
neous evolution, we make two types of fits: 

"Regge" : The 'Regge factorisation' approach discussed in Sect. 2, where 
there is no direct Pomeron contribution and no inhomogeneous 
term in the evolution equation. 

"pQCD" : The 'perturbative QCD' approach discussed in Sect. 3, where 
these effects are included. 

We make separate fits to the recent HI LRG (prel.) [1] and ZEUS Mx [3] 

CT?^"''' data, applying cuts > 3 GeV^ and Mx > 2 GeV, and allowing 
for overall normalisation factors of 1.10 and 1.43 to account for proton 
dissociation up to masses of 1.6 GeV and 2.3 GeV respectively. Statistical 
and systematic experimental errors are added in quadrature. The strong 
coupHng is set via as{Mz) = 0.1190. We take the input forms of the 
DPDFs at a scale Q§ = 3 GeV^ to be 

zYP{x^,z,QI) = frixr) C^z^-{\ - zf\ (10) 
zg^[xy,z,Ql) = /p(xp) CgZ^-{\ - zf^, (11) 

where /p(.xp) is given by (4), and where ap(0), Ca, Aa, and Ba (a ~ q,g) 
are free parameters. The secondary Reggeon contribution to the HI data 
is treated in a similar way as in the HI 2002 fit [1], using the GRV pi- 
cnic parton distributions [17]. Good fits are obtained in all cases, with 
xVd-o.f. = 0.75, 0.71, 0.76, and 0.84 for the "Regge" fit to HI data, 
"pQCD" fit to HI data, "Regge" fit to ZEUS Mx data, and "pQCD" fit to 
ZEUS Mx data respectively. The "pQCD" fits are shown in Fig. 4, includ- 
ing a breakdown of the different contributions. The DPDFs are shown in 
Fig. 5. Note that the "pQCD" DPDFs are smaller than the corresponding 
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Figure 4. "pQCD" fits to (a) HI LRG and (b) ZEUS Mx data. 



"Regge" DPDFs at large z due to the inclusion of the higher-twist 7£P — > qq 
contribution. Also note that the "pQCD" DPDFs have slightly more rapid 
evolution than the "Regge" DPDFs due to the extra inhomogeneous term in 
the evolution equation (9). There is a large difference between the DPDFs 
obtained from the HI LRG and ZEUS Mx data due to the different 
dependence of these data sets; see also [7,8]. 

The predictions from the two "pQCD" fits for the charm contribution 
to the diffractive structure function as measured by ZEUS using the LRG 
method [18] are shown in Fig. 6. Our HI LRG fit gives a good description, 
while our ZEUS Mx fit is too small at low /3. Note that the direct Pomcron 
contribution is significant at moderate (5. These charm data points were 
included in the determination of DPDFs from ZEUS LPS data [2], but only 
the resolved Pomeron i^^* (f — *■ cc) contribution was included and not the 
direct Pomcron (7*P cc) contribution. Therefore, the diffractive gluon 
distribution from the ZEUS LPS fit [2] needed to be artificially large to fit 
the charm data at moderate (3. 
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Figure 6. Predictions for ZEUS LRG difFractive charm production data using DPDFs 
from the "pQCD" fits to (a) HI LRG and (b) ZEUS Mx data. Note the large direct 
Pomeron (7*P — > cc) contribution at moderate /3. 



5. Conclusions and outlook 



To summarise, diffractivc DIS is more complicated to analyse than inclu- 
sive DIS. CoUinear factorisation holds, but we need to account for the direct 
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Pomeron coupling, leading to an inhomogeneous evolution equation (9).'^ 
Therefore, the treatment of DPDFs has more in common with photon PDFs 
than with proton PDFs. The HI LRG and ZEUS Mx data have a differ- 
ent dependence, leading to different DPDFs. This issue needs further 
attention. For a NLO analysis of DDIS data, the direct Pomeron terms, 
C2,p, and Pomeron-to-parton splitting functions, Pav, need to be calculated 
at NLO. There are indications [16] that there are large 7r^-enhanced vir- 
tual loop corrections ('K-factors') similar to those found in the Drell-Yan 
process. As with all PDF determinations, the sensitivity to the form of 
the input parameterisation, (10) and (11), and input scale Qg needs to be 
studied. The inclusion of jet and heavy quark DDIS data, and possibly 
^D(3) measured [22], would help to constrain the DPDFs further. 

The extraction of DPDFs from HERA data will provide an important input 
for predictions of diffractive processes at the LHC. 
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